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Abstract 

^^ By applying density functional theory, we find strong evidence for an itinerant nature of magnetism in two families of iron pnictides. 
'- ' Furthermore, by employing dynamical mean field theory with continuous time quantum Monte Carlo as an impurity solver, we 
Z^ observe that the antiferromagnetic metal with small magnetic moment naturally arises out of coupling between unfrustrated and 
frustrated bands. Our results point to a possible scenario for magnetism in iron pnictides where magnetism originates from a strong 
TT^Jnstability at the momentum vector (n, n, n) while it is reduced by quantum fluctuations due to the coupling between weakly and 
^^ strongly frustrated bands. 
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Since the discovery of high-Tc superconductivity in 
LaOFeAs 1 1], great effort has been devoted to pursuing higher 
superconducting transition temperatures Tc L2J as well as to the 
'understanding of the various phase transitions in these materi- 
als as a function of temperature, pressure and doping Ij, |4|, |51] . 
Density functional theory (DFT) calculations |6, 7] show that 
electron-phonon coupling - though non-negligible - is not 
strong enough to explain the high Tc observed in these systems. 
Instead, magnetically mediated pairing has been proposed to 
account for the superconducting state due to its proximity to a 
stripe-type antiferromagnetic (AF) metallic state ||8i|9[]. 

However, the origin of the stripe-type AF metal, i.e., whether 
■it comes from the itinerant nature of the Fermi surface (FS) 
"nesting r9l ll0llllLll2lll3n or a localized picture of exchange in- 
teractions between local spins 11141. 15 , 16] , is still under debate. 
The scenario of FS nesting was severely challenged by recent 
DFT calculations [17], which found a disconnection between 
Fe moment and FS nesting. However, a recent revision of the 
nature of magnetism in the iron pnictides - also within DFT, but 
with more precise optimized structures - has reestablished the 
close connection between itinerancy and magnetism 1131 . In 
this work, we will present further clear evidence of the itinerant 
nature of magnetism by performing DFT calculations for a few 
families of iron pnictides. 

On the other hand, the disagreement llisl 11911 of the mag- 
nitude of the Fe moment in the stripe-type AF metal between 
experimental observations and DFT calculations based on ex- 
perimental structures is still controversely discussed. Vari- 
ous mechanisms have been proposed from DFT calculations. 
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for example, negative effective on-site electronic interaction 
U |20l,l21], a low moment solution within GGA+U stabilized 
by the formation of magnetic multipoles |22, 23] or frustration 
between local spins 1 14 1 15112412511 . Recently, also mean-field 
calculations based on a five-band Hubbard model with positive 
U obtained a small magnetic moment comparable to experi- 
mental results i26ll27ll28ll . However, the existence of various 
sets of DFT-derived hopping parameters i29ll30ll31ll casts doubt 
on the reliability of the model parameters used in the model cal- 
culations. 

From DFT calculations it is well known ISl [iSl ll9|] that the 
physical properties of the iron pnictides are highly sensitive to 
details of the structure as well as to details of the exchange and 
correlation (XC) functionals. For instance, structural optimiza- 
tion within the local spin density approximation (LSDA) leads 
to almost perfect agreement |18] with the most recent experi- 
mental value of the Fe moment II32I1 . This indicates, that quan- 
tum fluctuations, which are only insufficiently incorporated in 
the common approximations to DFT, could strongly improve 
the agreement with experiment. 

To further explore this, we employ dynamical mean-field 
theory (DMFT) 13311 combined with continuous time quantum 
Monte Carlo (CTQMC) Q simulations. A feature fliat is com- 
mon to different sets of DFT-derived hopping parameters is the 
presence of some strongly and some weakly frustrated Fe Jid 
bands ll29[l30ll3lll . In the present work we consider a minimal 
two-band Hubbard model which should capture this feature by 
considering one band with frustration and the second one with- 
out frustration and we investigate the effect of the coupling of 
these two bands on the magnetism of the system. We would like 
to remark that while a realistic description of the Fe systems 
needs at least a five-band model, the present two-band model 
should be sufficient for the proposed analysis. We find that an 
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AF metallic state is present in a wide range of the interaction 
parameter U when one band is highly frustrated and the second 
one unfrustrated, while the state is absent when both bands are 
equally frustrated [35]. 

Our results from DFT and DMFT suggest that a strong insta- 
bility at a momentum vector (tt, n, n) is the promising mech- 
anism for the itinerant magnetism observed in iron pnictides 
while quantum fluctuations originating from the coupling be- 
tween weakly and strongly frustrated bands reduce the mag- 
netic moment and make it more comparable to the experimental 
observations. 



fulfills the rotational invariance condition of L'^ = U' +2J and 
ignore the spin-flip and pair-hopping processes. The operators 
are written in the standard notation of the multi-band Hubbard 
model 11351]. In order to solve this model we employ the two- 
sublattice DMFT method ll33ll which includes the local quan- 
tum fluctuation effects and can account for the AF state, com- 
bined with CTQMC simulations L34J. Existing DMFT studies 
combining with DFT calculations are focused on paramagnetic 
state ||43[ 144 l45l l46l 14711 . Our calculations are performed on the 
Bethe lattice. 



2. Method and Model 



3. Results and Discussion 



In order to quantify the FS nesting and hence the itinerant 
nature of the magnetism, we calculate (i) the q-dependent Pauli 
susceptibility at w=0 with the constant matrix element approx- 
imation, defined as 
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where a and jS are band indexes and q and k are momentum vec- 
tors in the Brillouin zone, and (ii) the k^ dispersion of the FS. 
These calculations were performed with the full potential lin- 
earized augmented plane wave method as implemented in the 
WIEN2k code Q with RK^^^ = 7. While 40000 k points ai-e 
used in calculating the k~ dispersion of the FS, a three dimen- 
sional grid of 128 X 128 X 128 k and q points are employed for 
the susceptibility. All calculations were performed in the scalar 
relativistic approximation. 

In order to perform the DFT analysis, (i) we use the avail- 
able experimental structures for LaOFeAs ll37ll . BaFe2As2 11381] 
and (ii) we obtain fully optimized structures for BaFe2-,iCoj:As2 
within the virtual crystal approximation and for a few hy- 
pothetical compounds Uke LaOFeSb and BaFe2Sb2 using the 
Car-Parrinello [39] projector-augmented wave [40] method. 
Our optimized structures compare well with the experimental 
ones il3ll4lll . Part of our results are double-checked by the full 



potential local orbital (FPLO) method f42]. Results are consis- 
tent among these methods. Throughout the paper, the Perdew- 
Burke-Ernzerhof generalized gradient approximation (GGA) to 
DFT has been used. 

For our model calculations we consider the following two- 
band Hubbard model 
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where f,„ (f^,) is the intra-band hopping integral between 
nearest-neighbor (next nearest-neighbor) sites with band in- 
dices m - 1,2. For simplification, we neglect inter-band 
hybridizations. U, U' and J, are the intra-band, inter-band 
Coulomb interaction and Ising-type Hund's coupling, respec- 
tively. In our calculations we set U' 
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Figure 1 : (Color online) Comparison of normalized q-dependent Pauli suscep- 
tibilities calculated within GGA at fixed q- = n along the [110] direction be- 
tween arsenide and antimonide of (a) 1111 compounds and (b) 122 compounds. 
The normalization factors are the susceptibilities of the coiresponding arsenide 
systems for each type of compound at qo = (0, 0, ti). Please note that the peak 
position is not exactly at q;^ = (n, n, n) since the electron and hole FSs are nearly 
nested rather than perfectly nested. 



Table 1: Magnetic moment calculated within spin-polarized GGA based on 
experimental structures of LaOFeAs and BaFe2As2 and optimized structures 
of the hypothetical compounds LaOFeSb and BaFe2Sb2. 



m(fiB) 



LaOFeAs LaOFeSb 



1. 



2.2 



BaFe2As2 BaFe2Sb2 
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2.5 



In Fig. [T] we present the comparison of q-dependent Pauli 
susceptibilities between arsenide and antimonide 1111 and 
122 compounds normalized by the susceptibilities of the cor- 
responding arsenide systems for each type of compound at 
qo - (0, 0, 7t). Here we only consider q, - n since in both 
LaOFeAs and BaFe2As2 the Fe spins are AF ordered along c as 
experimentally observed and we show the results for q^ - qy. 
We observe that in all these four compounds a strong peak 
around (\„ - (n, n, n) is found, which supports the presence of 
stripe-type magnetically ordered states as observed experimen- 
tally. Most importantly, in contrast to an earlier DFT study il7|] 
where the magnetic moment increases as As is replaced by Sb 
in both LaOFeAs and BaFe2As2 while the susceptibilities at q^r 
decrease in BaFe2Sb2 compared to BaFe2As2, which is inter- 
preted to be evidence of disconnection between FS nesting and 
magnetism, our results show that in both compounds, the sus- 
ceptibilities at ({„ increase as the replacement takes place and 
simultaneously the magnetic moment is also enhanced as dis- 
played in Table. [T] This discrepancy is attributed to the im- 
proved precision of optimized lattice structures in our DFT cal- 



culation 11311 . Therefore, the close connection between itiner- 
ancy and magnetism remains. 
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Figure 2: (Color online) (a) Normalized q-dependent Pauli susceptibilities of 
LaOFeAs and BaFe2As2 at fixed g, = q^ = n along the [001] direction cal- 
culated within GGA. The normalization factors are the susceptibilities of each 
compound at q^ = (n, n, 0), and the experimental structure is used, (b) Normal- 
ized q-dependent Pauli susceptibilities of BaFe2-.vCojAs2 at fixed qi= tt along 
the [110] direction calculated within GGA. Here .v = 0,0.2,0.5. The normal- 
ization factors are the susceptibilities of each compound at qo = (0, 0, tt). 

In Fig.|2](a), we show the q-dependent Pauli susceptibilities 
of LaOFeAs and BaFe2As2 at fixed q^ - n and qy - n, normal- 
ized by the susceptibilities of each compound at q,', - (n, n, 0). 
We find that while a tiny increase of the susceptibilities is found 
from q,') to q„ in LaOFeAs, indicating a dispersionless FS along 
c and therefore nearly perfect two dimensional physical prop- 
erties, a stronger enhancement is seen in BaFe2As2 suggesting 
a three dimensional FS topology as shown in Fig. [3] (a), (b), 
even though it is a layered compound. Such a three dimension- 
ality of the FS has been proposed to be the mechanism for a 
nearly isotropic critical field in (Ba,K)Fe2As2 [48]. The com- 
mon feature in the susceptibilities for LaOFeAs and BaFe2As2 
is that the values at q„ are larger than those at q^ which can 
account for the AF arrangement of Fe spins along c, although 
the interlayer interaction is believed to be small. Such a con- 
sistency again implies a close relation between itinerancy and 
magnetism. 
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Figure 3: (Color online) k- dispersion of FS calculated within GGA as a func- 
tion of Co-doping in BaFe2-.vCovAs2 around F at .v = (a), at Jt = 0.2 (c), at 
X = 0.5 (e) and around X at .r = (b), at x = 0.2 (d), at jc = 0.5 (f). 

It has already been shown that changes in the FS topology 
can explain the different nature of structural and magnetic phase 
transitions in the 122 compounds underpressure illll49tl . In the 
following we will further demonstrate that the phase transitions 
in BaFe2As2 under Co-doping are also related to the change 
of FS topology. It is known from experiments 15011 that for 
BaFe2-.vCo^-As2 at x = 0, the system shows a stripe-type AF 



metal. At x = 0.2 the system becomes superconducting, and 
finally at x = 0.5, both magnetization and superconductivity 
disappear. In Fig. [3] we show the kr dispersion of the FS for 
these three cases around F (left panels) and X (right panels). 
We find that the FS along c becomes more dispersive with Co 
doping. At X = a strong instability in the Pauli susceptibility - 
though the FS nesting is not perfect - is present at q„ as shown 
in Fig.|2](b) solid curve, supporting the AF state. At x = 0.2, 
the distortion of the FS indicates a strong suppression of the 
q„ instability (see Fig.|2](b) dashed curve) and therefore of the 
magnetization. However, the superconducting state which, ac- 
cording to the spin fluctuation theory |51], is related to the in- 
stabilities around q„ may in such a situation be more favourable 
than magnetization. At x = 0.5, the FS around F shrinks to a 
Fermi pocket and is even more distorted around X, suggesting 
that no obvious instabihty in the susceptibility will be present 
(see Fig. |2](b) dotted curve) and leading to the disappearance 
of both magnetic ordering and superconductivity. 

After having presented various evidence for the itinerant na- 
ture of magnetism in iron pnictides, we will investigate in 
what follows a possible mechanism for the reduced magnetic 
moment observed experimentally compared to DFT calcula- 
tions ill[ll8lll9tl . As mentioned in Section [1] we would like 



to extract essential physics from a simplified model as intro- 
duced in Section |2] In this model, local quantum fluctuations 
are included in the calculations by employing DIMFT. Fig. |4] 
shows the magnetization as a function of U at two tempera- 
tures. Combining these results with the analysis of density of 
states in Ref. [35], we conclude that when two bands are equally 
highly frustrated (see Fig.|4](a)), a first-order phase transition 
from a paramagnetic metal to an AF insulator state is observed 
and an AF metallic state is absent, while, if one band is unfrus- 
trated and the second one highly frustrated (see Fig.|4](b)), sev- 
eral continuous phase transitions appear separately in these two 
bands, and an AF metal with small magnetic moment appears. 
This indicates that an AF metallic state with small magnetic 
moment emerges out of the coupling between highly frustrated 
and unfrustrated bands, which is the case in iron pnictides as 
mentioned in Section [T] rather than out of a pure frustration 
effect. 
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Figure 4: (Color online) Magnetization per band of a two-band Hubbard model 
at two temperatures calculated by DMFT(CTQMC) as a function of interaction 
strength U with (a) two bands equally highly frustrated and (b) one band highly 
frustrated, the other unfrustrated. 



4. Conclusion 

In summary, we presented various evidence for the close re- 
lation between itinerancy and the magnetism observed experi- 
mentally. Our results support the itinerant nature of magnetism 
in iron pnictides and suggest that a strong instability at the nest- 
ing vector q„ is responsible for the magnetism in iron pnictides. 
We also propose that the reason why the reduced magnetic mo- 
ment in the iron pnictides cannot be reproduced by LSDA or 
GGA calculations employing the experimental lattice structure 
is the insufficient incorporation of quantum fluctuations in such 
calculations. By applying the DMFT approach to a simpli- 
fied two-orbital Hubbard model, we propose that coupling of 
strongly frustrated and unfrustrated bands may be the mecha- 
nism for the reduction of the magnetic moment, rather than a 
pure frustration eff'ect. 
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